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ABSTRACT 


Fly ash, a waste material of the thermal power plant is being 
produced enormously in our country, whereas its utilization is 
very limited. At present it is generally disposed off near the 
vicinity of thermal power plant in the form of lagoon ash. In the 
present investigation liquefaction analysis of embankment 
constructed by " Upstream method of stage construction" is carried 
by empirical methods, pseudo static method, steady state method 
and simplified dynamic analysis. 

In this investigation fly ash samples at Satpura thermal 
power plant (STPP) Sarni (M. P.) and Hansdeo Thermal Power Plant 
(HDTPP) Korba west (M.P.) sites were tested to establish necessary 
geotechnical properties required for the evaluation of 
liquefaction potential of fly ash disposal sites. As an example, a : 
typical staged ash embankment consisting of five stages (total | 
height of 80 ft. (24.3m)) has been evaluated. ; 

(1) Empirical analysis was performed based on the earthquake i 
severity Index (ESI) procedure as outlined by ICOLD applicable for : 
earthquake zone of the site and typical upstream embankment i 

I 

configuration. | 

(2) Pseudo static analysis was carried out as per the ; 
Recommendation given in 1.3.1893:1984. 

(3) In steady state method , consolidated undrained test with pore | 
pressure measurement were performed on samples compacted at 85 %, j 

j 

90% , 95 % and 100 % Standard Procter density and on a very | 

i' 

loose sample and the void ratio versus Sus relationship were j 
established. Settled ash slurry tube samples were prepared in the ' 
laboratory at various consolidation pressure simulating the field 
condition and a relationship between void ratio versus 



consolidation pressure was established. Octahedral stresses likely 
to be developed in different layers of settled ash, is computed 
from the applied vertical consolidation pressure on the soil 
surface. Corresponding to that pressure, the void ratio was 
determined from the tube sample results. The corresponding 
undrained steady state shear strength (Sus) values were determined 
from Sus versus void ratio relationship. Embankment stability 
analysis was cai'ried out using PCSTABL5 program and it was found 
that the embankment is safe for five stages of upstream 
construction . 

(4) Simplified dynamic analysis was carried out using SHAKE 
program. The shear wave velocity and shear modulus were determined 
from the standard correlation using the laboratory test results. 
The analysis gave a maximum vertical deformation of 1.2 
inch. (3.048 cm) for STPP and 1.07 inch. (2.71 cm) for HDTPP 
embankment site whereas the empirical relationship gave the 
vertical deformation in the range of 1.92 inch (4. 88 cm.) to 0.58 


inch (1.47 cm) . 
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CHAPTER 1 

Introduction 


In India demand of electrical energy is increasing day by 
day. Due to ample coal reserve in India, it has resulted in the 
construction of coal fired power plants. Most of the power plants 
being constructed or have been proposed are coal fired power 
plants. As a result the production of power plant waste fly ash is 
increasing day by day. Fly ash is a by product of those thermal 
power stations which are using pulverized coal as a fuel. At 
present 65 thermal power plants are operating in India, producing 
about 40 million tons of ash as a waste material in our country 
per annum. With the commissioning of the super power plants and 
increasing use of low grade coal of high ash content, the present 
ash production is likely to be tripled. This will pose a serious 
disposal and ecological problem in addition to occupying large 
tracts of scarce of cultivatable land. 

The beneficial usage of ash in concrete , brick making ,soil 
stabilization treatment and other applications have been 
recognized; only a small fraction of total ash produced in India 
is currently utilized in such applications . Until full utilization 
of ash is made, ash produced from the power plant will be disposed 
off on ground. Since the volume of ash to be disposed off is 
large, technical skill is essential for its disposal. It can be 
disposed off either by dry scheme or by wet scheme. In India ash 
from the most of the power plants is disposed off by wet scheme in 
the vicinity of the power plant as a waste material . In wet 
disposal scheme fly ash is generally transported as a slurry 
through pipe and disposed off in an impoundment (starter dyke). 
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When the starter dyke is approximately full four basic options are 
available for the disposal of ash: 

1) Constructing new lagoons using either conventional 

construction material 

2) Having the fly ash in existing lagoons hauled to another 
disposal site 

3) Raising the existing dyke using conventional material 

4) Raising the dyke using the fly ash excavated from the lagoon. 
The fourth option is considered to be very cost effective because 

fly ash used for constructing dyke would in addition to saving the 
earth filling cost, enhance the disposal capacity of the 
lagoon. This can be achieved by any of the following method : 

1) Upstream method of stage construction 

2) Downstream method of stage construction 

3) Centreline method of stage construction 

Upstream method : 

In this method an initial starter dyke is constructed at the 
downstream toe of proposed embankment and then embankment 
construction progresses upstream in stages as the height of the 
embankment increases. This technique requires minimum volume of 
structural fill for the raising dykes and it has better 
environmental pollution control compared to other methods because 
the constructed face becomes the final face of the constructed 
embankment. The major drawback of this method is that a part of 
embankment is constructed on the uncompacted ash, which is not 
structurally sound material. Probability of liquefaction of 
hydraulically filled lagoon ash is very high when it is subjected 



Downstream Method of Stage Construction 



CentreLine Method of Stage , Construction 


Fig. 1.1 : Different Methods of Stage Construction 



Downstream method : 


In this method, construction is carried out on the downstream 
side of the starter embankment no that the crent of tlu^ dam in 
shifted progressively towards the downstream. The main advantage 
of this method is that none of the embankment is constructed on 
the uncompacted ash as in case of the upstream method. Major 
drawback of this technique is that it requires relatively large 
volume of fill to raise the embankment. 

Centreline method : 

It is basically a type of downstream method in which crest of 
the embankment is not shifted in the downstream direction but 
raised vertically upward above the crest of the starter dam. This 
method has many design construction, environmental and operational 
problems and as such it is generally not adopted. 

The present thesis covers the seismic design aspect of new 
tailings dam or fly ash embankments constructed by upstream method 
of staged construction. 

In chapter 2 embankment stability under seismic condition has 
been reviewed for both compacted and uncompacted embankments. More 
emphasis has been given on compacted embankment constructed over 
the settled ash slurry. 

In chapter 3 experimental setup and its result has been 
discussed. 

In chapter 4 theoretical analysis has been carried out. 
Chapter 5 presents results and discussion. Conclusion and 
recommendation for future work has been given in chapter 6 . 



CHAPTER 2 

Literature Review 

2.1 General : 

Liquefaction is a phenomenon wliere in a saturated sand and 
other granular materials such as fly ash looses a large percentage 
of shear resistance due to monotonic or cyclic loading and flows 
like liquid until the shear stress acting on the mass are as low 
as its reduced shear resistance. Earthquake causes strong 
reciprocating, horizontal (most damaging) and vertical forces due to 
inertia of the structure such as earth dam, tailing dams etc. , 
and the ground motions during earthquake are periodic but highly 
irregular. The cyclic stress reversal between compression and 
extension as well as repeated compressive loading tends to 
decrease continuously the volume of the particulate material which 
gives rise to increase in pore pressure .Hence the effective 
strength of the saturated sand decreases. The liquefaction 
potential and degree of liquefaction will depend on type and 
relative density of the sand and the average total stress as well 
as the level and duration of cyclic and repeated loading. 

2.2 Materials Susceptible to liquefaction: 

When an element of saturated soil is subjected to shaking 
(stimulated in laboratory by cyclic load for earthquake) the pore 
pressure increases .The pore pressure generated into the soil due 
to cyclic stress depends on a numbers of factors, but on the most 
elementary level, soil type can be considered the fundamental 
parametre .For saturated cohesive soils subjected to undrained 
cyclic loading, pore pressure is generated but usually to a 
limited degree. As a result the static undrained shear strength 
of a clay may be reduced somewhat by cyclic loading; but the soil 
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retains its basic integrity and liquefaction does not occur. 
Consequently clayey noils are not UMually ol. primary concern in 
determining embankment stability under earthquake loading 
conditions . 

Another class of material is represented by very coarse 
permeable deposits of cohesionless material. High permeability of 
such materials prevents undrained conditions from developing 
during cyclic shear; in effect pore pressure generated by cyclic 
loading dissipate as fast as they are generated. Here too 
liquefaction susceptibility is seldom of major concern. 

Dense sand still represents another class of material .Most 
commonly such materials are in centreline or downstream type 
embankments constructed of fly ash placed in lifts and compacted 
under controlled conditions . Pore pressures are generated in 
response to undrained cyclic loading and may momentarily reach 
high values. However the tendency of such material to dilate 
quickly which reduces the pore pressure, causing the material to 
stabilize. Even after a large number of cyclic reversals the 
cumulative pore pressure is limited, although the sample may 
experience some degree of permanent deformation. Severe loss of 
strength does not occur under these conditions and therefore the 
compacted material does not usually experience severe seismic 
instability. 

The most problematic class of materials consist of loose 
saturated sands of relative density of 30-50%. This category may 
apply to most upstream embankments constructed on uncompacted fly 
ash deposits .When it is subjected to strong and/or extended 
seismic shaking the undrained cyclic loading conditions is likely 
to occur. Pore pressures are generated during each applications of 
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cyclic stress reversal, and the effects become the cumulative 
during the duration of shaking. Under the application of uniform 
cycles of cyclic stress pore pressure may reach the Icv-ol of (:lu' 
applied confining stress. At this point the effective stress 
within the foundation material becomes momentarily zero and 
initial liquefaction is said to have occurred. Following initial 
liquefaction, the embankment may rapidly undergo excessive strain 
only with a few additional cycles of stress. 

Seed (1979) had reviewed the factors related to the 
nature of material and its depositional characteristics governs 
the liquefaction susceptibility. Many of those factors suggest 
that the liquefaction susceptibility of tailings may be greater 
than -that of natural soil deposits at the same relative density. 
Those factors include grain characteristics and gradation, method 
of deposition, aging, previous strain history and over 
consolidation ratio. These factors suggest the comparison between 
natural soil deposits and tailings deposits at equivalent relative 
densities must be made with caution and that tailings samples for 
laboratory testing must be prepared in a way that reflects their 
depositional and environmental character in situ. With regards to 
failure of tailings dam due to ground rupture is remote and it may 
be possible only when the dyke is built over the actual fault 
zone. This type of failure is generally reviewed at the time of 
site selection and it shall not be, therefore, addressed here. By 
far the common problem in upstream method of tailings dam des.ign 
is to ensure that the dyke will be stable under the anticipated 
levels of seismic shaking. The purpose of the literature review is 
to provide an overview of procedures for assessing fly ash 
tailings dyke stability under earthquake loading. 
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The technology of seismic analysis of earth dam stability is 
a young one . Only within last three decades the mechanism has been 
developed for an understanding of actual soil behaviour during 
earthquakes and those mechanism to conventional earth dam is not 
simple. In case of tailings dam, their unique material properties 
further complicate the problem. For a complete seismic analysis of 
dam /tailings dam due to earthquake loading requires {l)the 
estimation of seismic acceleration, (2) the response of 
soil/material to cyclic loading and (3) seismic analysis. 

2.3 Analysis of compacted embankment : 

Dams constructed of clay on clay or rock foundations have 
withstood extremely strong shaking under acceleration of .35 -.5 g 
with no permanent damage (Seed et.al. 1977, 1978). Rockfill dams 
have also performed well. Dams of this class are unlikely to 
experience failure due to seismic liquefaction. In the context of 
tailing embankments, they include downstream and water retention 
type embankments constructed of clay, rockfill, mine waste or of 
other material that has been well compacted under controlled 
conditions. For embankment of this class, seismic analysis are 
commonly performed using either a pseudo static or deformation 
approach . 

2, 3. 1 Pseudo-static Method ; 

In this analysis the problem is treated as a static one, with 
a horizontal force applied by earthquake expressed as the product 
of seismic coefficient multiplied by tlie weight of the poteuLial 
sliding mass . 

Pseudo static methods have long been used and they remain 
most useful for seismic stability analysis in ‘cases where cyclic 
liquefaction or major pore pressure build up is not anticipated. 



Nevertheless, the pseudo static approach does not accurately model’ 
true embankment behaviour under major earthquake shaking. As such 
pseudo- static methods must be considered largely empirical in 
nature. However, this method is routinely used in practice for 
compacted embankments . 

2.3.2 Deformation approach : 

Deformation approach may be applied to embankment materials 
that do not experience severe strength reduction during 
application of cyclic stresses, in effect this method assumes that 
the embankment slope will not undergo complete failure by sliding 
during an earthquake but that the embankment will undergo some 
permanent deformation. The basis for evaluating the embankment 
behaviour under seismic shaking is whether or not the computed 
deformations lie within acceptable limits i.e. to prevent 
overtopping. 

Procedures originally developed by Newmark{1965) have ben 
used to estimate embankment deformation during earthquake shaking 
.This procedures use limiting equilibrium methods to estimate 
permanent displacement experienced by a rigid-body failure mass, 
other simplified techniques proposed by Makadisi and Seed (1977) 
include consideration of the dynamic response of the embankment . 


2.4 ANALYSIS OF UNCOMPACTED EMBANKMENTS ; 

Embankments of uncompacted material, usually unconsolidated 
settled fly ash, pose an entirely different analytical problem 
than the compacted one. In uncompacted material the pore pressure 
can be expected to result from strong seismic shaking, and the 
material can no longer be assumed to retain all or the most of its 



original static strength. Upstream embankments of hydraulically- 
discharged tniling.s nrn the moat problematic variety of 
uncompacted embankments , but tlii.c., cralegcty a 1 mo Inc.-ludoM 
downstream and centreline type embankments that are saturated 
within significant portions and are constructed of uncompacted 
cycloned tailings . 

Many types of techniques, ranging from very crude to 
extremely complex, are available for seismic stability analysis of 
uncompacted embankments, but most suffer from their theoretical 
limitations or practical drawbacks. In a hierarchy of complexity, 
these methods include empirical evaluations, pseudo static method, 
in conventional or modified form, steady state method, simplified 
liquefaction analyses, and complete dynamic analysis. 

2.4.1 EMPIRICAL APPROACH : 

On the most simplistic level, it is instructive to make a 
preliminary assessment of seismic stability on the basis of the 
observed performance of hydraulic fill embankment during 
earthquakes .Catastrophic seismic failures or near failures of 
uncompacted hydraulic fill embankments have been evaluated in the 
context of the level of earthquake shaking they have experienced 
-The review of the literature suggests that most upstream 
embankments have survived acceleration up to about 0.15g (where g 
is the acceleration due to gravity) . 

2. 4. 1.1 Performance of Embankment Dams : 

Behaviour of embankment dams under seismic condition include 
settlement, horizontal movement, cracking, pore pressure build up, 
slope slumping and failure, internal erosion, increase of seepage, 
and reservoir and impoundment breaching in the extreme cases . 

The Seismic crack movement of embankment dams reflects dam 
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deformation associated with compression, lateral spreading and 
slope movements induced by earthquakes. Bureau et. al . (1985) 
reviewed seismic performance of earth and rock fill dams including 
the crest movement (settlement and horizontal movement) . They 
introduced an Earthquake Severity Index (ESI), which reflects the 
influence of both the intensity and duration of the earthquake and 
is defined as : 

ESI = A(M-4.5)^ Eq. 2.1 
where A = Peak ground acceleration at the dam site ; 

M = Earthquake magnitude 

Further an empirical relationship is given between the ESI and 
the crest settlement expressed as a relative settlement ratio 
(crest settlement divided by height) .This correlation plot updated 
with data from Lo and Klohn is shown in Figure 2.1. a The solid 
line represents the interpolated trend through all the data 
points . 

Figure 2.1.b presents a similar plot, also updated by data 
from Lo and Klohn (1992) correlating the horizontal relative 
crest movement (horizontal crest movement divided by dam height) 
with the Earthquake Severity Index, ESI. 

2.4,2 PSEUDO STATIC METHOD : 

2.4. 2.1 Conventional Pseudo Static Method : 

This method is also used in practice for uncompacted tailings 
embankments but they have even less theoretical validity than in 
their conventional application to compacted embankments because 
pore pressure due to cyclic loading are not accounted for. Failure 
might have been correctly predicted by using a seismic coefficient 
corresponding to the actual maximum acceleration experienced by the 
embankment during the earthquake- But there is uncertainity in 
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selection of appropriate strength values, hence to make a 
successful prediction is very rare. In general , this method can be 
adopted for those cases where seismic risk is so low that it can be 
reasonably certain that significant pore pressures due to cyclic 
loading will not develop. As a precaution total stress parametre 
may be used for saturated uncoinpacted material for low to moderate 
permeability in order to account in acrude way for pore pressures 
due to shear generated by the rapidly applied eartquake loading. 
2.4.2. 2. Modified Pseudo Static Method : 

Several investigators have attempted to address defiencies in 
the conventional pseudo state approach by accounting for pore 
pressures induced by cyclic shear in a more explicit way. One 
approach is to first consolidate the sample in the laboratory and 
then subject it to a number of cycles of undrained shear 
corresponding to the magnitude and duration of the design 
earthquake, and finally to shear the sample in a static undrained 
mode (assuming that liquefaction has not already occured) .A 
strength envalope results that is similar to a conventional 
envalope from static consolidateed tests , except that the pore 
pressures induced by the cyclic shear are incorporated. Lee and 
Roth (1977) has described the use of such a strength envalope in 
the analysis of hydraulic fill dam. 

2.4,3. STEADY STATE METHOD ; . 

The steady state of deformation for any soil mass (Poulos et . 
al. 1985) is that state in which the mass is continuously deforming 
at constant volume , constant normal effective stress , constant 
shear stress and constant rate of shear strain. The steady state 
of deformation can be achieved only after the structure is 
completely remoulded and all particle orientation effects have 



reached a steady state condition and after all particle breakage, 
if any is complete , so that the shear stress and the rate of 
deformation remain constant. The steady state can be reached in 
drained or uyndrained shear. 

2.4. 3.1 Steady state line 

Steady state line of a soil is most commonly expressed by a 
best fit line plotted between void ratio vs. effective minor 
principal stress (on log scale) or deviator stress. 

2. 4. 3. 2. Liquefaction evaluation procedure : 

Liquefaction potential in steady state method is determined 
by stability analysis which requires the determination of shear 
strength and the shear stress in situ .The shear strength used in 
the liquefaction analysis is the in situ undrained steady state 
shear strength which is difficult to determine accurately because 
it is very sensitive to void ratio. 

Poulos et al . (1985) has recommended the detailed step wise 

procedure for liquefaction evaluation. The steps in the analysis 
related to the measurement of in situ strength are listed below -. 

1) Determination of in situ void ratio 

2) Determination of steady state void ratio or density as 
a function of effective stress using compacted specimen. 

3) Determination of undrained/ steady-state strengths for 
undisturbed specimen. 

4) Correction of measured undrained steady state strengths to in 
situ void ratio. 

5) Calculation of in situ driving shear stress and the factor of 
safety against liquefaction. 

2.4. 3. 3 Determination of in-situ void ratio ; 

For clays, it is relatively simple to measure in situ void 
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ratio using the procedure developed by Hvorslev (1949) , but for sand 
or fly ash it is difficult to get a correct in situ void ratio. 
Generally fixed piston sampling, frooTiing of Iho ground and 
coring or block sampling in the test pit can provide satisfactory 
result but each one requires extra care and are not practicable 
for deeper strata. 

It is very difficult to measure in situ void ratio 
accurately at depth in a bore hole , particularly for a cohesion 
less soil. The in situ undrained steady state strength is highly 
sensitive to void ratio, therefore the in situ void ratio has to 
be measured with a procedure like that described in subsequent 
steps to obtain a suitable measurement of the undrained steady 
state strength (Poulos et al . 1985). 

2. 4.3.4. Determination of steady state void ratio as a function 
of effective stress using compacted sample : 

A procedure is required for correcting the laboratory 
measured steady state strengths to the in situ void ratio because 
the undisturbed samples of loose sands have always lower void 
ratio in laboratory than in situ. Even if a sample is obtained by 
freezing the ground and coring, the density increases 
significantly during reconsolidation in the laboratory to in situ 
stresses . Sometimes the increase in density during consolidation 
is so high that it causes soil to become dilative in laboratory 
whereas it was contractive in situ. This means that a sample gets 
converted from a state in which it is in liquefiable to a state in 
which it is not liquefiable. 

The influence of apparently minor change in void ratio on 
undisturbed steady state strength is shown in Figure 2.2. It shows 
the results of measurements of steady state strength for a clean 
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sub angular sand, quartz and tailings sand as a function of 
density. As shown in this figure even a slight change in the void 
ratio significantly effects the Sus value and a correction 
procedure has been suggested for obtaining in situ void ratio. 

The procedure for correcting laboratory measured undrained 
steady state strength to the in situ void ratio is based on two 
observations: (1) the slope of steady state line on the semi log 
plot is affected chiefly by the shapes of the grains in a given 
soil and (2) the vertical position of the steady state line is 
affected by even small differences in grain size distribution 
(Castro 1975) . 

The correction procedure currently used requires that the 
steady state line be obtained by testing five or six compacted 
specimens of identical soil. The soil used should be the same as 
that in the ground and at this stage grain shape distribution 
needs be well preserved, since it has the major effect on the 
slope of the steady state line. The results of the tests on 
compacted specimens are plotted on the "state diagram" (Figure 
2.3) a plot of void ratio versus effective minor principal stress 
a-^'.The best fit line through the points representing the steady 
state is the steady state line (SSL) . Since it is the steady state 
shear strength that is needed for liquefaction analysis, it is 
convenient to plot the results of the undrained triaxial tests in 
terms of void ratio versus undrained steady state shear strength 
on the failure plane, Sus, as shown in Figure 2.4.b rather than in 
terms of cr as shown in Figure 2.4 a. The Sus from the results of 

3s 

each consolidated undrained triaxial is computed test using the 
following equation : 
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Sus = q cos(j)s 


Sin(ps 


O'-, + q 

3k 


( cr - Aus ) + q 

3r 


3s 


Eq . 2.2 
Eq . 2.3 

Eq . 2.4 


in which = Principal stress difference 

d- 3 s = Effective minor principal stress at the steady state 
= Effective minor principal stress at start of shear 

Aus= Pore pressure induced in the test specimen at the steady 
state of deformation. 

4>s = Steady state friction angle 

The quantities q and Ausare obtained directly 

during triaxial tests. 

The steady state line obtained from the compacted samples 
is used to correct the strengths of "undisturbed " samples , 
determined in the following steps. 

2.4. 3.5 Determination of undrained steady state strengths for 
undisturbed samples : 

A series of consolidated undrained triaxial test is performed 
on "undisturbed " samples to determine the average steady state 
strength. It is best to insure that the sample is contractive in 
nature. One procedure that can be used is to consolidate the 
undisturbed specimens to high effective stress, but not so high 
that the correction needed in section 2. 4. 3. 6 becomes excessive. 
Steady state can still be determined at lower effective stress but 
the accuracy of the measurement is poorer than the contact ive 
samples . Soils having some plasticity , other tests must be used to 
determine the steady state line, such as rotation shear test and 


vane shear test . 
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The steady state point (point s in Fig. 2.5 ) f or each 
undisturbed sample is plotted on the same steady state diagram , as 
shown in Fig. 2. 4. a or 2.4.b with the steady state line for 
compacted specimen. The vertical distance between the individual 
point and the steady state line for compacted specimen is assumed 
to be due chiefly to minor differences in grain size 
distribution. The undrained steady state strengths shown in 
Fig. 2. 3. were obtained at the void ratio after consolidation ,not 
at the in situ void ratio .therefore correction of the results to 
the in situ void ratio must be made as described in the next 
steps . 

2. 4. 3. 6 Correction of measured undrained steady state strengths 
to in situ void ratio : 

In situ void ratio is measured during taking out the 
"undisturbed" samples .Using the void ratio of undisturbed samples 
measured in laboratory after consolidation, it is joined with the 
in situ void ratio keeping this line parallel to the steady state 
line of the compacted samples on the state diagram. The estimated 
in situ undrained steady state strength is selected from the 
abscissa, as shown in Fig. 2. 6 . 

2.4. 3.7 Determination of in situ driving shear stress and the 
factor of safety : 

The in situ driving shear stress in the zone being evaluated 
is calculated by conventional methods of stability analysis.lt is 
the minimum shear stress that is required to maintain the static 
equilibrium. The shear stress along the potential failure surface 
in the zone being evaluated can be computed using limiting 
equilibrium approaches, such as method of the slices, circular arc 
type , wedge type or by using a non circular surface .One can make 
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(After Poulos et. al. 1985) 
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Fig. 2.6 Correction of Measured Undrained Steady State Strengih 
FOR Difference between In Situ Void Ratio During Test 



type , wedge type or by using a non circular surface .One can make 
an assumption about the distribution of stresses along the failure 
surfaces to calculate this shear stress. 

Using the stability analysis , the driving shear stress "x " 

d 

is computed. If its average value is less than the undisturbed 
steady state shear strength in all zones along trial failure 
surface , then liquefaction cannot occur. 

The factor of safety against liquefaction , is given by 

Undrained steady state shear strength 

F = 

^ Shear stress required to maintsin static equelibrium 

The undrained steady state strength calculated above is the 
minimum strength that a contractive soil mass can have at the void 
ratio in situ. Therefore the factor of safety need only be 
sufficiently greater than 1 to ensure that any uncertainty in the 
void ratio is accounted for. 

2.4.4 Simplified Dynamic analysis;; 

2. 4. 4.1 Estimating Dam And Embankment Earthquake Induced 
Deformation : 

A simple yet rational approach to the design of dam and 
embankment under earthquake loading is based on the concept of 
permanent deformation as proposed by Nevmiark(1965) but modified 
to allow for the dynamic response of the dam and embankment as 
proposed by Seed and Martin (1966) rather than the rigid body 
behaviour.lt assumes that the failure occurs on a well defined 
slip surface and that the material behaves elastically at stress 
levels below failure but develops a perfectly elastic behaviour 
above yield. The method involves the following steps : 

1) Determination of a yield acceleration at which a potential 



sliding' surface would de'velop a factor of safety of unity is 
determined . 

2) Determination of earthguake induced acceleration in the 
embankment using the dynamic response analysis. 3) For a given 
potential sliding mass, when the induced acceleration exceeds 
the calculated yield acceleration, movements are assumed to 
occur along the direction of the failure plane and the 

magnitude of the displacement is evaluated by a simple double 

integration procedure. 

2. 4. 4. 2 Determination of yield acceleration : 

The yield acceleration (ky) is defined as the average 

acceleration producing a horizontal inertia force on a potential 
sliding mass so as to produce a factor of safety of unity and thus 
cause it to experience permanent displacement . This can be 
achieved using one of the available methods of stability analysis. 
In the present study the modified Bishop method has been used to 
calculate the yield acceleration for a circular slip surfaces 
using a pseudo static analysis. 

After determining the yield acceleration for a certain 

location of the slip surface, the next step in the analysis is to 
determine time history of earthquake induced average acceleration 
for the particular sliding mass. 

2.4.4. 3 Determination of earthquake induced acceleration: 

In order to determine the permanent deformation for a 
particular slip surface, the time hi-story of onrthqunko Ludiu'od 
average acceleration must first determined. As shown in Fig. 2. 7 at 
each time step the forces acting along the boundary of the sliding 
mass are calculated from the orresponding normal and shear stress 
of the finite elements along that boundary . The resultant of these 
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n 

F(t)= Et ( t ) L, + Ohi ( t ) d| 
i=I '"'I 

n = number of elements along the sliding surface. 

ka^(t) = F(t)/w 

Figure 2.7 Calculation of Average Acceleration from 
Finite Element Response Analysis 



forces divided by the weight of the sliding mass would give the 
average acceleration , kav (t) , acting on the sliding mass would 
give the average acceleration. Comparing the time history of crest 
acceleration with that of average acceleration for different depth 
of potential sliding mass, the similarity in the frequency content 
is readily apparent (it generally represents the first natural 
period of the embankment) while the amplitudes are shown to 
decrease as the depth of the sliding mass increases towards the 
base of the embankment. The maximum crest acceleration is 
designated by Umax and kmax is the maximum average acceleration 
for a potential sliding mass extending to a specific depth y. 
Fig. 2. 8 shows a relationship showing the variation of the maximum 
acceleration ratio kmax/umax with depth for a range of embankments 
and earthquake loading conditions.lt would then be sufficient for, 
design purposes to estimate the maximum crest acceleration in a 
given embankment due to specified earthquake and use this 
relationship to determine the maximum average acceleration for any 
depth of the potential sliding mass. 

Shake (Earthquake response analysis of horizontally layered 
site of version 1988) program has been used in this study to 
determine the variation of maximum acceleration witli depth using 
some published results of response computations. 

Fig 2 . 8 shows the results of variation of maximum average 
acceleration ratio with depth based on the results reported by 
Ambraseys and Sarma (1967) and Seed et.all (1978 on the basis of 
finite element calculation) .Knowing the yield acceleration and the 
time history of average acceleration for a potential sliding mass, 
the permanent displacement can be calculated using the fig. 2.9, 
showing the variation of y/h with ky/kmax.for different magnitude 






Fig. 2.8 Variation of Maximum Acceleration 
Ratio With Depth of Sliding Mass 



Fig. 2.9 Variation of Yield Acceleration With (a) 
Permanent Displacement (b) Average Normalized Displacement 


Normalized 




of earthquake . 

2,4.5 COMPLETE DYNAMIC ANALYSIS ; 
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Dynamic analysis represents the state of the art for seismic 
analysis of uncompacted embankments subjected to strong shaking. 
Dynamic embankment analysis for conventional dams are described by 
Seed et al . (1975) Lee and Roth (1977) and Makdisi et al . (1978) Finn 

and Byrne (1976) Doery and Byrne (1982 ) describe dynamic analysis 
procedures applied to centreline type tailings embankments . In a 
dynamic analysis of a centreline embankment Lo et al . (1982 ) discuss 
the effects of assumed slimes liquefaction on stability of the 
embankment . 

Seed(1979 b) and Finn ( 1982 } summarize procedures by performing 
dynamic analyses of conventional dams respectively. In general, 
dynamic analyses incorporate the following steps : 

1) Determine the preearthquake static stress in the embankment. 

2) Determine the dynamic response of the embankment and foundation 
to base rock motions. 

3) Evaluate the dynamic soil behaviour 

4) Evaluate embankment stability . 

2.4. 5.1 Determine the preearthquake static stress in the 
embankment : 

This is usually performed using a static finite element 
analysis to determine the initial effective normal stress and 
shear stress along the potential failure .These stresses are 
essential for proper interpretation of laboratory cyclic data. 
2.4,5. 2 ; Determine the dynamic response of the embankment and 
foundation to base rock motion ; 

On the basis of a well instrumented earthquake in a similar 
geologic environment and of similar size and duration , a bed rock 



acceleration time history is established. Based on this base 
motion time history , the response of the embankment is determined 
by dynamic finite element modeling , using either equivalent linear 
or non linear procedures. The end result is to determine the 
acceleration time history at every node in the finite element 
mesh. These time histories must be converted to an equivalent 
series of uniform cyclic stress application. 

2 . 4. 5. 3 ; Evaluate the dynamic soil behaviour : 

This is ordinarily accomplished by performing cyclic stress 
tests for conditions corresponding to those indicated by the 
results of static stress analysis with the laboratory results 
corrected to account for several factors that influence the in 
situ behaviour .Failure in laboratory samples may be defined in 
terms of initial liquefaction or with reference to a particular 
level of axial strain. 

2. 4. 5. 4 : Evaluate embankment stability : 

Stability of the embankment may be addressed in terms of 
either stresses or strains occurring locally within the 
embankment. The stability of the embankment is evaluated by 
comparing stresses induced by the earthquake to those required to 
cause liquefaction or to exceed prescribed liinits of local strain. 
The locations within the embankments for which induced stresses 
exceed those required to cause liquefaction or produce 
unacceptable strain determine the extent of liquefaction or 
unacceptable performance. Embankment deformations can be assessed 
qualitatively on the basis of strain potential for the individual 
elements, which corresponds to the strain to that would be 
experienced if the element were not constrained by surrounding- 
soil. Elements with excessive strain potential may be considered 



to offer no excessive resistance to slide movement within the 
embankment. Under these conditions, a stability analysis can be 
performed to indicate the overall factor of safety against sliding 
through zones of liquefaction or excessive deformation (Seed et.al. 
1975) . Alternatively, cyclic induced pore pressures for individual 
elements, together with effective stress strength parametres, can 
be incorporated into a limit equilibrium analysis to compute a 
post earthquake factor of safety. In addition, Lee and Roth 
(1977) recommended to calculate total embankment deformation on the 
basis of gravity loads and softened material properties to 
determine whether they lie within the acceptable limits. 

2.5 Packages used : 

Two software packages have been used in this investigation, 
they are 

1 . PCSTABL5 

2 . SHAKE 
2.5.1 PCSTABL5 ; 

it is a personal computer based program developed by Siegel 
(1975) and Boutrup (1977) in FORTRAN IV language for the general 
solution of slope stability problem s by two dimensional limit 
equilibrium methods, the factor of safety against instability of 
a slope is performed by the method of slices , modified Janbu 
method and Spencer method (1967) for slip surfaces of any shape. 

The PCSTABL5 features a random .search for gonornrion of 
potential failure surfaces for subsequent determination of more 
critical surface and their factor of safety. One technique 
generates circular surfaces , another surfaces of sliding 
character and third more general irregular surfaces of random 



shape. This program also handles the heterogeneous soil system, 
anisotropy soil strength properties, excessive pore pressure due 
to shear, static ground water and surface water , pseudo static 
earthquake loading and surcharge boundary loading. 

2. 5. 2 SHAKE : 

Shake, a personal computer based program is developed by Seed 
H. B. Schnabel B. and Lysmer J. (1972) which is later modified by 
Sun J . and Lai S. It computes the responses in a system of 
homogeneous, viscoelastic layers of infinite horizontal extent 
subjected to vertically travelling shear waves. The program is 
based on the continuous solution to the wave equation (Kanai 
1951) adapted for use with transient motions through the Fast 
Fourier Transform algorithm (Cooley and Tukey, 19650. The non 
linearity of the shear modulus and damping is accounted for by the 
use of equivalent linear soil properties (Idiris and Seed 1968 
Seed and Idris 1970) using an iterative procedure to obtain values 
for modulus and damping compatible with the effective strain in 
each layer. The following assumptions are made in this analysis : 

1) Soil system is infinitely extended in the horizontal direction 

2) Each layer is completely defined by its value of shear modulus, 
critical damping ratio density and thickness . 

3) The response in the system are caused by the upward propagation 
of shear waves from the underlying rock formation. 

4) The shear waves are given as acceleration values of equally 
spaced time intervals . Cyclic repetition of the acceleration time 
history is implied in the solution. 

5) The strain dependence of the modulus and damping is accounted 
for by an equivalent linear procedure based on an average 
effective strain level computed for each layer. 



The program is able to handle systems with variation in both 
moduli and damping and takes into account the effect of the 
elastic base. The motion used as a basis for the analysis , the 
object motion, can be given in any one layer in the system and new 
motions can be computed in any other layer. 

layered site using different earthquake record data such as 
Pasadena earthquake (1952) , Kerntaft earthquake etc. program has 
been used in this study to determine the variation of maximum 
acceleration with depth using some published results of response 
computations. Following sets of operation can be performed by 
this program : 

1) By reading the input motion it determines the maximum 
acceleration, scales the values and computes the predominant 
period . 

2) By reading the soil deposit properties , it computes the 
fundamental period of the deposit. It computes the maximum stress 
and strains in the middle of each sub layer and modifies the value 
of modulus and damping with a specified percentage of maximum 
strain . 

3) It computes the new motion at the top of any sublayer inside 
the system or outcropping from the system which can be printed, 
plotted or punched as per the requirement . . 

3) It also computes the response spectra for the motion. 

4) Amplification function between any two sublayers can be 
computed, printed and plotted. 

;i 

5) Without changing the predominant period or duration of the 
record time interval can be increased or decreased. 

6) Computed motion can be set as a new motion. 

7) Stress or strain time history in. the middle of any sublayer can 



be computed, printed and plotted. 

These operation can be performed by exercising the 
available options in the program. 


various 



CHAPTER 3 


Experimental Details 

3.1 General : 

This chapter includes the source of material, details of 
their properties , preparation of the sample, type of test 
performed, test details and test results.. 

3.2 Source of Material : 

Three types of flyash samples, obtained from three different 
sites, have been investigated in this program. Samples were 
obtained from Satpura Thermal Power Plant (STPP) Sarni dist. Betul, 
M.P., Korba Super Thermal power Plant (KSTPP) Korba dist. 
Bilaspur; M.P. and Hansdeo Thermal Power Plant (HDTPP) Korba west 
dist. Bilaspur M.P.. All the samples were obtained from ash 
disposal ponds of those sites. 

3.3 Properties Of Fly ash ; 

The physical and chemical behaviour of fly ash is 
governed by several factors such as coal source, degree of coal 
pulverization, design of boiler unit, loading and firing 
conditions and handling and storage methods. A change in any of the 
above factors can result in detectable changes in the ash produced 
and it's properties 

Thus it is not surprising that a high degree of variability 
can occur in ash properties not only between power plants but 
within a single power plant. The degree to which any change 
affects the utilization potential of ash is a function of nature 
and degree of changes, and the particular application for which 
the ash might be utilized. In this investigation following 



distribution, specific gravity, optimum moisture content, maximum 
dry density, shear strength parameters, and steady state 
strength . 

Grain size distribution curve: 

3 ^ ^®P^®s®^ts the combined results of the sieve analysis 
and hygrometer analysis. The grain size distribution analysis 
gives an idea about type of material and gradation. It is clear 
from the Fig. 3.1 that the samples used are s.ilt to sand sized 
particles . 

Specific gravity: -■ 

Specific gravity tests were done in the laboratory with the 
help of specific gravity bottle. Specific gravity of KSTPP fly ash 
reported by Was an. (1995) has been used here. The results are 
given in the Table 3.1. 

Standard Procter Test: - 

Standard Procter test was carried out to determine the 
maximum dry density and optimum moisture content of STPP and HDTPP 
ash samples. Standard Procter test result of KSTPP sample given by 
Wasan (1995) has been used here, these results are given in Table 

3.1. 


Shear strength Parameter :- 

For saturated ash samples, effective 'strength parameters were 
determined by consolidated undrained triaxial test with pore 
pressure (CU) measurement. STPP samples have been tested at 90% 
and 95% maximum dry density, whereas HDTPP sample has been tested 
at 95% maximum dry density. For KSTPP the shear strength parameter 
given by Wasan (1995) is used in the investigation. The shear 
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at 95% maximum dry density. For KSTPP the shear strength parametre 
given by Wasan M. (1995) is used in the investigation. The shear 
strength parametres (cohesion (c' ) and effective friction angle 
(<!>')) for all the three sites are summerized in the Table 3.1. 



Table 3 . 1 Geotechnical properties of samples 

Sample 1. STPP Fly ash 
Specific Gravity 1.98 

Optimum Moisture Content 33 % 


Maximum Dry Density 

1.11 gm/cc 

Shear 

strength Parametre 


C' 

} 

1 (95% M.D.D. ) 

0 

<t>' 

35.5° 

C 

> 

0 

\ 

<p' 

1 (90% M.D.D. ) 

33.7° 


Sample 2. HDTPP Fly Ash 


Specific Gravity 


Optimum Moisture Content 


1.94 


Maximum Dry Density 

1.145 gm/cc j 

Shear 

strength Parametre 


C' 

> 

0 

<l>' 

1 (95% M.D.D. ) 

31.5° 


Sample 3 . KSTPP Fly ash 
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3,4 Sus Test : 

In order to determine the steady state strength of various 
ash consolidated undrained triaxial test with pore pressure measureme 
has been carried out on samples of 38.1 mm. diametre and 76.2 mm. 
high sample. The back pressure equal to cell pressure of 1.5 
Kg/cm is used to saturate the sample .To ensure the complete 
saturation Bishop s parametre ( B = — ) was checked and B 

3 

nearly equal to unity ensured practical limit of 100% saturation. 
After saturation, the sample has been consolidated for eight hours 
at cell pressure equal to 2.5 Kg/cm^and back pressure of 1.5Kg/cm^ 
i. e. effective consolidation pressure of 1 kg/cra^and then the 
sample was sheared using constant rate of shear strain of 6 
mm/min. . In case of non cohesive soil such fly ash is no pore 
pressure lag is likely to occur at the applied strain rate of 
6mm/min..The deviatoric stress is calculated after applying the 
area correction from the proving ring reading and the strain in 
the sample is measured from a dial gauge. It has been observed 
that the ash sample attains its steady state generally at 15-16 % 
strain. At steady state the shear stress, deviator stress and 
effective minor principal stress become more or less constant with 
further the increase of strain. Graphs were plotted between effective 
principal stress versus strain and deviator stress versus strain 
and value of q and O' were taken from these graphs to calculate 
the S 

US 

3.4.1 Test Sample : 

Three types of samples have been prepared for consolidated 
undrained test with pore pressure measurement. They are as 
follows : - 

1) Compacted sample 

2) Tube sample 



3) Loose sample 


Compacted sample; 

Based on standard Procter test results, samples were prepared 
at 85%, 90%, 95% and 100 % Standard Procter dry density in 
miniature Procter compaction sampler. Fly ash of required weight 
has been taken in china bowl and mixed thoroughly with required, 
quantity of water . Then it is compacted in miniature Procter 
sampler in three layers and each layer was compacted twenty five 
times. For STPP no test has been carried out at 100% maximum dry 
density because it has been observed that samples tested is 
dilating at 90% and 95% of Standard Procter dry density shows a 
dilating tendency and therefore it was decided to not to test it 
as it will not likely to liquefy. 

Tube sample: 

Fly ah was thoroughly mixed with tap water ( 1:10) to make it 
in slurry form and then laboratory consolidated fly ash samples 
were prepared by sedimentation. The slurry was poured in a tube 
fitted with porous stone covered with filter paper at top and 
bottom and allowed to consolidate under its own weight. The excess 
water is allowed to escape taken out through the bottom porous 
stone and after that a porous stone fitted with filter pater at 
top and bottom is placed over the sample and a dsired load is 
placed over it through a hanger arrangement Sample is then allowed 
to consolidate for next few days. Samples so formed was 
transferred into the miniature compaction sampler carefully to 
allow minimum disturbance and was cut in standard length by means 
of sharp edge. The STPP samples were consolidated in the tube at 
0.3 kg/cm^, 0.6 kg/cm^and 0.9 kg/cm^ pressure and HDTPP samples 
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were prepared at consolidation pressure of. 6 kg/cm^, 0.9 kg/cm^, 
and 1.2 kg/ctn^ . 

Loose sample: 

Undisturbed samples taken from field generally get disturbed 
during sampling trasporting and testing. Also it is very 
difficult to prepare compacted sample at less than 85% MDD or at 
low consolidation pressure. In order to get good quality of loose 
sample for testing, it was decided to prepare loose sample in the 
cell itself, using vaccuum technique. 

3.5 Test Results : 

Sust ests were performed on the Procter compacted samples, 
tube consolidated samples and the loose sample. Fig. 3. 2, Fig. 3. 3 
and Fig. 3. 4 show the variation of deviator stress versus strain of 
compacted samples of STEP, HDTPP and KSTPP samples, respectively. 
From these figures it is clear that deviator stress increases with 
dry density of sample. The deviator stress also increases with the 
strain and eventually it becomes constant at a very high strain. 

Fig. 3. 5 to Fig. 3.7 shows the variation of change in pore 
pressure with starin for compacted samples of STEP, KSTPP and 
HDTPP, respectively. It is clear from those results that the 
variation in pore pressure with strain is much higher in loose 
samples than in dense samples and the variation of pore pressure 
gradually decreases as the density increases. Also at higher 
strain, change in pore pressure becomes more or less constant. In 
case of STPP sample change in pore pressure for 90% and 95% sample 
{Fig. 3. 5) is found to be negative at higher strain i.e. samples 
are dilating at those densities. 

Fig. 3. 8 to Fig. 3.10 show the variation of effective minor 
principle stress versus strain for the compacted samples of STPP, 



Deviator stress (Kg/cm^2) 
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strain 


Fig. 3.2 variation of Deviator stress with strain (STPP) 
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Fig. 3.6 variation of Pore Pressure with strain (KSTPP) 




Change in Pore Pressure (Kg/cm^2) 



Fig. 3.7 variation of pore pressi 



Effec 



strain 


Fig. 3.8 Variation of [EMPS] with strain (STPP) 
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KSTPP and. HDTPP samples. Theses results are mirror image of the 
pore pressure variation versus strain, off course the magnitude is 
different . 

Fig. 3.11 and Fig. 3. 12 show the variation of deviatoric 
stress versus strain of STPP and HDTPP tube samples respectively. 
It is clear from these figures that the deviator stress increases 
with increaes of consolidation pressure. Fig. 3.13 and Fig. 3. 14 
show the variation of change in pore pressure with strain of STPP 
and HDTPP samples respectively. Fig. 3. 15 and Fig. 3. 16 show the 
variation of effective minor principal stress versus strain of 
STPP and HDTPP tube samples. These test resultsa figure indicate 
that the behaviour of tube samples is dilative in nature. Low 
confining pressure during the test is probably one of the reasons 
for the dilative nature of the sample. 

Steady state shear strength has been calculated from 
the (CU) test results. 

SSL For Ash sample : 

The Sus values for all the test results are calculated from 


the equation 2.2 taking the difference of principal stress, 
effective minor principal stress and change of pore pressure at 
steady state. Figs. 3.2 to Fig. 3.4 gives the difference of 
principal stress at steady state for compacted and loose samples 
and figs. 3.11 to Fig. 3. 12 for tube samples. The change in pore 
pressure is taken from Fig. 3.5 to Fig. 3.7 for STPP KSTPP and 


HDTPPsample, respectively and cfrcctive mi not; pi i nci pnl .Ml i cu?;? al 
steady state is taken from Fig. 3. 8 to Fig 3.10for compacted and 
loose ash samples. Similarly for tube samples the respective 
values were taken from Figs. 3. 13 to 3.16. In Table 3.2 Sus value 


and corresponding void ratio of compacted samples and loose 
samples are summerized. Sus value of tube samples are summerized 



.6 Kg/cm^2 
.9 Kg/cm^2 
1 .2 Kg/cnn^2 



Fig. 3.1 1 Variation of Deviator Stress vs. Strain (HDTPP) 
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Fig. 3.16 Variation of [EMPS] with strain (STPP) 



Table 3.3 Sns and void ratio for tube consolidted fly ash samples 

Cons o lidation 0 h q Water Void Sun Effective 

Pressure .Content Minor 

Kg. /Cm ' 1 Degree I Kg/Cm 'I After Ratio, Kg/Cm^' Principal 

Test , Stress 

Kg/cm^ 


For STPP Fly Ash 


32.1 

2. 97 

54 . 1 

1 . 1 


34 . 0 

3.52 

52.6 

1.08 

2.92 

34.5 

3.86 

51.0 1 

1.05 

3 .18 


For IIDTP Fly Ash 


33.2 3.29 


53.2 1.055 


48.7 


47.3 




























SSL FOR CPMPACTED SAMPLE (STPP) 
SSL FOR CPMPACTED SAMPLE (HDTPP) 


SSL FOR CPMPACTED SAMPLE (KSTPP) 
O’ SSL FOR TUBE SAMPLE (STPP) 

♦ SSL FOR TUBE SAMPLE (HDTPP) 




CHAPTER 4 

Theoretical Analysis 

4.1 General: 

This chapter deals with the determination of Sus, shear 
modulus, and shear wave velocity at different depth. 

4.2 Determination of steady state strength at different depth : 

Sus value at field void ratio can be determined by correcting 
the Sus value of the undisturbed samples tested in the laboratory 
as explained in section 2. 4. 3. 6 .Due to non availability of field 
samples here an attempt has been made to study the liquefaction 
phenomenon of tailings dam by interpolating results between the 
compacted dense samples and loose samples prepared in triaxial 
cell itself. Tube sample results as shown in Fig. 3. 13 & 3.15 of 
STPP and Fig. 3.14 & 3.16 of HDTPP clearly indicate the dilative 
behaviour and hence there is no possibility of liquefaction to 
occur in such types of samples. In the present analysis, 
therefore only the results of compacted and loose samples have 
been used. 

As discussed in section 2.4.3 Sus value is function void 
ratio. 

As the field void ratio at different depth is not presently 
available, it has been decided to get the corresponding field void 
ratio from the void ratio vs. consolidation pressure results of 
tube samples . As shown in Fig .4.1 a straight line relationship 
is found between void ratio and log of octahedral consolidation 
pressure. For obtaining field void ratio first octahedral pressure 




Fig. 4.1 Variation of octahedral consolidation pressure 
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is determined at different depth using the correlation between 


effective stress and poission ratio as given in Eq. 4.1 
ll±2y)(r' 2 . 

cr = :r = -TT-Cr' 

0 3 3 

cr^ = Octahedral pressure 


Eq. 4.1 


cr' = Effective Pressure 
V = Poission Ratio = .5 (Assumed) 

From the octahedral pressure calculated at different 


depth, 

corresponding void ratio is determined by using Fig. 4.1. Using 
this void ratio the corresponding Sus value has been determined 
from SSL diagram of compacted and loose samples at different 
depths, from the Fig. 3. 17. The estimated Sus values at different 
depth is tabulated in Table 4.1 to 4.2 respectively, for STPP and 
HDTPP settled ash layers. 

4.3 Determination of shear modulus shear wave velocity at different 
depth : 


Shear modulus of soil at any depth can be calculated by using 
the formula given by Hardin-Black (1968) , later modified by 
Christiano (1980) , which is as follows : 


Gmax 

Where 


17.760 (2.973 
(1 + e) 



) 


0.5 


G 

2.7 



Gmax = Shear Modulus (KSF) 
e = Void ratio 

cr^ = Octahedral pressure (PSF) 

G = Specific Gravity 


Shear wave velocity at any depth can be calculated by using 
the following formula : 

{ Gma X -yO . S 
V = ( p 1 



where : 


V = Shear Wave Velocity 

Gmax = Shear Modulus 
p = Mass Density 

Void ratio as detrmined in section 4.2 is used to determine 
the shear modulus and shear wave velocity at different depth. 
Using the above formula shear modulus and shear wave velocity 
different depth has been calculated and the results are given in 
the Table 4 . 1 & Table 4.2 for STEP and HDTPP, respectively. 



Tabls 4.1.. Calculation of Sus Sfisax* Modul'^s and. sh.eaxr ' 
velocity at different depth of embankment made of STPP Fly Ash 


DEPTH 

EFFECTIVE 

STRESS 

OCTAHED- 

RAL 

PRESSURE 

VOID 

RATIO 

SHEAR 

MODULUS ■ 

SHEAR 

WAVE VEL. 

Sus 

m 

cr' 

Kg/ cm^ 

a 

0 

Kg/cm^ 

e 

6max ' 

Kg/cm^ 

V 1 

m/Sec) 

Kg/cm^ ; 

*0.760 

0.106 

0.070 

0.920 

102.063 

85.042 

1.67 

3.048 

0.337 

0.225 

0.920 

182.332 

194.194 

1.67 

5.334 

0.409 

0.273 

1.093 

154.641 

186.404 

0.293 

6.096 

0 .440 

0.293 

1.090 

161.110 

190.263 

0.299 

7.620 

0.502 

0.334 

1.084 

173 .486 

197.436 

0.336 

9.144 

0.563 

0.375 

1.079 

185.220 

204.004 

0.352 

10 .668 

0.624 

0.416 

1.075 

196 .410 

210.075 

0.368 

12.192 

0.686 

0.457 

1.071 

207.131 

215.733 

0.395 

13.716 

0.747 

0.498 

1.067 

217.442 

221.038 

0.413 

15.240 

0.808 

0.539 

1.063 

227.391 

226.038 

0.423 

16.760 

0.870 

0.580 

1.060 

237.018 

230.773 

0.443 

18.288 

0.931 ^ 

0.621 

1,057 

226.352 

235.273 

0.453 

19.812 

0.993 

0 . 662 

1.055 

255.423 

239.565 

0.464 

21.336 

1.054 

0.703 

1.052 

264.254 

243.671 

0.486 ' 

22.860 

1.115 

0.744 

1.049 

272.863 

247.609 

0.508 


Conversion Factor 

1 Kg/Cm^ = 100 KN/m^ 







Table 4 . 2 . Calculation of Sus .Shear Modulus and shear 
velocity at different depth embankment made of HDTPP Fly Ash 


DEPTH 

EFFECTIVE 

STRESS 

OCTAHED- 

RAL 

PRESSURE 

VOID 

SlATIO 

SHEAR 

MODULUS 

SHEAR 

WAVE tEL. 

Sus 

m 

cr' 

Kg/cm^ 

cr 

Kg/cm^ 

e 

Gmax 

Kg/cm^ 

V 

Kg/ cm^ 

0.762 

0.109 

0.073 

0.900 

102.531 

83.923 

0.560 

3 .048 

0.364 

0.242 

0.900 

187.380 

203.067 

0.560 

5.334 

0.440 

0.293 

1.122 

147.055 

181.775 

0.165 

6.096 

0.468 

0.312 

1.108 

154.929 

186.578 

0.176 

7.620 

0.523 

0.349 

1.083 

170 .447 

195.699 

0.212 

9.144 

0.578 

0.386 

1.060 

185.691 

204.263 

0.249 

10.668 

0.634 

0.422 

1.039 

200.690 

212.352 

0.279 

12.192 

0.689 

0.459 

1.020 

215.469 

220.032 

0.299 

13 .716 

0.744 

0.496 

1.002 

230.048 

227.354 

0.313 

15.240 

0.800 

0.533 

0.986 

244 .443 

234.360 

0.352 

16.764 

0.855 

0.570 

0.970 

258.670 

241.083 

0.368 

18.288 

0.911 

0.607 

0.956 

272.741 

247.554 

0.404 

19.812 

0.966 

0.644 

0.942 

286.666 

253.795 

0.433 

21.334 

1.021 

0.681 

0.930 

300 .457 

259.823 

0.453 

22.860 

1.077 

0.718 

0.918 

314.120 

265.670 

0.474 


Conversion Factor 

1 Kg/Cm^ = 100 KN/m^ 


wave 



CHAPTER 5 


Results And Discussion 

5.1 General; 

This chapter deals with liquefaction analysis of a typical 
ash embankments constructed by upstream method as shown in Fig. 
5.1. For analysis geotechnical properties of STPP and HDTPP fly 
ashes, as discussed previously were used. The liquefaction 
analysis is carried out by empirical method, pseudo static method, 
steady state method, and simplified dynamic analysis following the 
procedures outlined in section 2.4 Analysis of Uncompacted 

Embankment " . 

5.2 Ash Embankment details 

Typically ash disposal dykes are constructed in several 
stages by "Upstream Method of Staged Construction" in five 
stages. As shown in Fig. 5.1 a typical upstream constructed 
embankment consist of five to six stages and consist of three to 
four types of soils i. e. foundation soil, compacted earthen 
starter dam and compacted staged construction of earth or ash 
material and settled ash. In this analysis for the simplicity it 
has been assumed that the geotechnical properties of foundation 
soil and compacted earth dyke is same. The other two materials are 
; compacted ash for stage II through stage V and settled lagoon 
ash. There are three types of soils assumed in this analysis are 
shown in Fig. 5.1. Soil I is foundation and starter embankment 
soil, soil II is the compacted ash used in the construction of 




Fig. 5.1 Upstream method of Staged Construction 



raising dykes and soil III is the settled Pond ash. The 
geotechnical properties used in the liquefaction analyses are 
summerized in the Table 5.1. These values are based on the test 
results and analysis as discussed in chapter 3 and chapter 4 . 


5.3 Liquefaction Analysis : 

5.3.1 Empirical approach 

As discussed in section 2. 4. 1.1, settlement of embankment 

depends on peak ground acceleration and earthquake magnitude . 
Different authors have given a correlation between peak ground 
acceleration and earthquake intensity as given in Table 5.2 . 

According to seismic zonal map of India (I . S . 1893 : 1984) both 
STPP and HDTPP sites fall somewhwere in the range of seismic zone 
II and III. The modified Mercelli intensities assciated with zone 
II and III are 6 and 7, respectively and also the earthquake 

magnitude for zone II and III are of the same order, i.e 6 to 7 . 
Earthquake magnitude of 6.5 is used here for the site under 

investigation . 

The value of peak ground acceleration is calculated from 
various empirical formulae given in Table 5 . 3 and the results are 

summerized- in the same Table . The ESI for each formula is 

calculated from the Eq. 2.1 and the percentage vertical settlement 
and horizontal movement are calculated from the relationship shown 
in Fig. 2.1. The absolute value of vertical settlement and 
horizontal movement for the configuration shown in Fig. 5.1 is 
also tabulated in the Table -. The range of vertical settlement is 
found to be .58 inch(1.46 cm.) to 1.92 inch(4.88 cm.)and the 

horizontal movement ranges from .25 inch(.62 cm.) to .91 inch(2.32 



Table 5.1 : Geotechnical Properties of soil and Fly Ash used to 

carry out liquefaction analysis 

Embankment Site HDTPP (KOREA WEST) 


Geotechnical 


Soil Type 


Properties 

I 

II 

III 

Specific 

Gravity 

2.65 

1.94 

1.94 

Moist 

Density 

1.71 

1.43 

0.84 (Dry) 

Saturated 

Density (gm/cc) 

1.79 

1.51 

1.36 

Shear strength 
Parametre 




Cohesion (C' ) 

0 

0 

0 

Ef f ctive 

Friction 

Angle {cp' ) 

30.0 

31.5 

26.0 

Shear Modulus 

1000 

Given in 

Table 4.2 

(KSF) 




Shear Wave 
Velocity 

800 

Given in 

Table 4 . 2 

Ft/sec 




Embankment Site 

STPP (SARNI) 


Geotechnical 


Soil Type 


Properties 

I 

II 

III 

Specific 

Gravity 

2.65 

1.98 

1.98 

Moist 

Density 

1.71 

1.39 

0 . 84 (Dry) 

Saturated 

Density (gm/cc) 

1.79 

1.50 

1.36 

Shear strength 
Parametre 




Cohesion (C') 

Eff ctive 

0 

0 

0 

Friction 

Angle {<()') 

30.0 

35.5 

26.0 

Shear Modulus 

lOOO 

Given in 

Table 4.1 

(KSF) 




-Shear Wave 




Velocity 

800 

Given in 

Table 4 . 1 

Ft/ sec 





1 





Table 5.3 - Deformation Calculation using Different Correlation as 
given in table 5 . 2 


Empirical 

Relation 

1 

ESI 

. 264 

MMI = 6 

AH 

c m 

1.47 

AD 

c m 

1.10 

ESI 

.552 

MMI = 7 

AH 

c m 

2.92 

AD 

c m 

1.24 

2 

. 520 

2.69 

1.27 

1.064 

4.88 

2.44 

3 

.496 

2.44 

1.24 

.992 

4.39 

2.18 

4 

.376 

1.96 

1.12 

.848 

4.24 

1.96 

5 

.432 

2.21 

1.17 

1.024 

4.60 

2.26 

6 

.456 

2.32 

1.19 

.816 

3.91 

1.91 
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cm. ) . It is to be noted that for the most conservative correlation 
by (Neumann 1954), the maximum settlement is found to be 1.92 
inch (4 .88 cm. ) . 

5.3.2 Pseudo Static analysis: 

Before carrying out the pseudo static analysis static 
stability analysis has been carried out to see whether it is 
statically safe or not. 

Tne seismic stability analysis is performed for pseudo static 
condition assuming suitable earthquake loading. The 
recommendations published by the Indian Standard (IS 1893 : 
1984) suggest that a seismic acceleration of . 05g acting in the 
horizontal direction may be used for the site. The factor of 
safety for static and dynamic conditionare shown in Fig. 5.1 
PCSTABL5; a software package has been used for both static and 
siesmic analysis which has been discussed in section 2.6. Using 
modified Bishop method, it has been found that the static factor 
of safety is 1.56 and 1.54 for STPP and HDTPP dyke site, 
respectively which is more than the required F.S. of 1.5 whereas 
under pseudo static condition factor of safety is 1.23 and 1.225 
for STPP and HDTPP dyke site, respectively. Various internatiionsl 
agencies have suggested that under pseudo static condition factor 
of safety should be of the order of 1 to 1.2. Computer print out 
for STPP dyke site has been given in the appendix. 

After ensuring that the dam is safe for both static and 
pseudo static conditions, further analysis has been carried out, 
assuming different values of horizontal acceleration to get a 
factor of safety of equal to one which indicates that the dam is 



unsafe for acceleration value higher than that at which factor of 
safety is one. It has been found that ash embankments both STPP 
and HDTPP are safe up to 0.11 g horizontal acceleration. 
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5.3.3 Steady State Method: 

In this method driving shear stress is determined using the 
undrained soil property of uncompacted ash. Driving shear stress 
is determined by trial and error using PCSTABL5' program for 
modified bishop method to obtan a factor of safety equal to one. 
Using a number of iteration it has been found that the driving 
stress for ash embankment at STPP site is 170 PSF and that of 
HDTPP site it is 175 PSF. From table 4.1 and 4.2, the minimum 
value of Sus is 600 PSF and 337 PSF for STPP and HDTPP which is 
much more than the driving shear stress . Hence there is no 

possibility of liquefaction to occur. The actual F.S.= 1.32 has 
been found for HDTPP dyke siteout from PCSTABL5 program for 

modified Bishop's failuer circles by substituting Sus values of 
various layers of 10 ft. (3.05 m increments). As it has already 
been discussed in the section 2.4.3 that Sus value is very 
sensistive to void ratio, this result is not very much 

satisfactory due to non availability of field void 

ratio. Nevertheless an attempt has been made to demonstrate, how 
the steady state method could be used for investigation of 

liquefaction potential of a fly ash dyke constructed by upstream 
method . 

5.3.4 Simplified Dynamic analysis: 

This analysis is another approach for assesing the 

liquefaction potential of a dam/tailing dams. In this method. 
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yield acceleration (Ky) , maximum average acceleration due to 
earthquake (Kmax) and time period (To) are obtained and from those 
results the maximum deformation of embankment is determined. 

As discussed in section 2. 4. 4. 2, yield acceleration (Ky) is 
the horizontal acceleration for which embankment has a factor of 
safety equal to one. This is calculated by pseudo static method as 
discussed in .the section 5.3.2. 

Maximum average acceleration and time period induced due to 
earthquake is determined using 'SHAKE' program for Pasadena and 
Kerntaft earthquake records. The reason for selecting the above 
earthquake records, is that the acceleration spectrum of these 
eatrhquake records are for medium dense soil . Similar one found 
at STTP and HDTPP sites .In order to calculate acceleration at 
different depth (to simulate earthquake condition) horizontal 
motion of 0 . 15g is given at the base of the embankment and 
induced acceleration at differnt depth is determined by using 
'SHAKE' program. It also gives the value of first natural period. 
As discussed in section 2. 4. 4. 3, the maximum average acceleration 
is determined by averaging the maximum acceleration from top of 
the embankment to the depth at which failure circle passes. The 
average of- these value gives the Kmax value. Knowing the value of 
Ky/Kmax and To (time period) maximum deformation of the embankment 
is calculated using the procedure outline in section 2.4.4. 
After Knowing the total deformation of the embankment one can be 
decide whether it is permissible or exceeding the limits. 

The shake analysis result is given in the appendix. The 
deformation analysis is carried out for three different height of 
the embankment, 60ft(18.29 ra) , 70 ft. (21.34m) and 80ft(24.38 m) 
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for both sites .The results are shown in the Table 5.3. As shown 
in the table for 80 ft. (24.38 m) , the maximum deformation from 
this method is found to be 1.07 inch (2. 72 cm) for HDTPP site and 
1.2 inch (3.05 cm) for STPP site. As mentioned in the section 
5.3.1, the deformation from empirical relation is found to vary 
between 1.92 inch (4.88 cm) inch to 0.58 inch(1.47 cm.) inch which 
is very comparable to those obtained for simplified dynamic 
analysis. One sample computer printout for HDTPP dyl<e site has 
been given in the appendix 



Table 5.4 Deformation calculation for embankment made by HDTPP 
Ash sample 


Embankment 

Time 

Period 

(Sec) 

Km a > 


Ky 

Ky /Km a x 

Um a X 

Um a X 

height 






Km ax gTo 

Cm. 

m 

K 

P 

K 

P 


K 

P 

K 

P 

K 

P 

24.38 

.94 

1.24 

. 188 

.186 

.11 

.585 

.59 

.013 

. 012 

2.26 

2.72 

21.34 

. 91 

1.22 

.165 

.170 

.11 

. 667 

.647 

. 006 

. 008 

0 . 89 

1.63 

18.29 

.86 

1.19 

.157 

.167 

.11 

.701 

.659 

.005 

. 007 

0.66 

1.37 


Deformation calculation for embankment made by STPP 
Ash sample 


Embankment 

Time 

Km a X 

Ky 


Um a X 


height 

Period 

(Sec) 








m 

K 

P 

K 

P 




K 

P 



24.38 

.98 

nra 

.183 

.185 

.11 

.601 

.595 

.014 

. 012 

2 . 11 

3 . 05 

21.34 

.94 

1.27 

.167 

.171 

.11 

.659 

.643 

.007 

. 007 

1.07 

1.78 

18.29 

.89 

1.22 

.157 

.164 

.11 

.701 

.671 

.005 

. 005 

0.69 

1.12 


K stands for Kerntaft earthquake rcord 
P stands for Pasadena earthquake record 
Umax - Settlement- 














































CHAPTER 6 

Conclusions 

In this investigation, liquefaction evaluation was carried 
out for two ash disposal sites satpura Thermal power plant (STPP) 
Sarni, and Hansdeo Thermal Power Plant (HDTPP) Korba West, 
located in Madhya pradesh. The essential geotechnical tests such 
as grain size distribution, specific gravity, standard procter 
test, triaxial test with pore pressure measurements and Sus test 
were performed which are required for the liquefaction analysis . 

The liquefaction analysis of a tailings dam is generally 
performed by empirical method, pseudo static method, steady state 
method, simplified liquefaction analysis and complete dynamic 
analysis. In this investigation, liquefaction analysis is carried 
out by first four method as stated above and the evaluation 
procedures and findings are summarized below: 

(1) Empirical analysis was performed based on the Earthquake 
Severity Index (ESI) procedure as outlined by ICOLD applicable for 
earthquake zone of the site and typical upstream embankment 
configuration. Range of vertical deformation by this method is 
found to vary between 1.92 inch. (4.876 cm) and 0.58 inch. 
(1.467cm) and the horizontal movement between .34 inch (.85 cm) to 
. 96 inch (2.44 cm) . 

(2) pseudo-static analysis was carried out as per the 
recommendation given in 1.8.1893:1984. Factor of safety equal to 
1.23 is found for horizontal seismic coefficient ofO.OSg which is 
recommended by I. S. 1893 :1984. 

(3) In steady state method, consolidated undrained test with pore 
pressure measurement were performed on samples compacted at 85 s , 
90 % ,95% and 100% standard Procter density and on a very loose 
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samples and the void ratio versus Sus relationship were 
established, settled ash slurry tube samples were prepared in the 
laboratory at various consolidation pressure simulating the field 
conditions and a relationship between void ratio versus 
consolidation pressure was established. The void ratio at 
different depth in a settled ash is estimated based on octahedral 
consolidation pressure likely to be developed in different layers 
of settled ash and corresponding to that pressure, void ratio 
were determined from the tube sample results. The corresponding 
Sus values were determined from Sus versus void ratio 
relationship. Liquefaction analysis was carried out and embankment 
is found to be safe against liquefaction for five stages of 
upstream construction. The factor of safety is found to be 
1.32 for HDTPP dyke site. 

(4) Simplified dynamic analysis was carried out using SHAKE 
program. Shear modulus and shear wave velocity were determined 
using the standard correlation using the laboratory test 
re suits. Maximum deformation of the five staged embankment 
constructed by upstream method of staged construction is 1.07 inch 
(2.72 cm) for HDTTP site and 1.2 inch (3.05 cm) for STPP site. 


Recommendation for Further Investigations: 

Due to the logistic problem and time limitations undisturbed 
sample could not be obtained from test sites. In this 
investigation entire exercise is carried out simulating field 
conditions in the laboratory, as discussed in the chapter two, the 
steady state undrained shear test strength and shear wave velocity 
are very susceptible to void ratio. It is therefore recommended 
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that similar investigations as presented here may be carried out 
on high quality of undisturbed samples obtained at various depth 
from ash lagoons. Complete dynamic analysis utilizing the finite 
element technique should also be persued with actual soil dynamic 
properties obtained from the ash lagoon sites. 



ll- V + t + 

S H A K E 5 

* * .♦ If: ♦ 4- 


A program for*. 

Earthquake Response Analysis 
or Horliontally Layered Sites 
Version of February 1988 

Originally authored by 
Schnabel, Lysmer, & Seed 
riodlfled by J. Sun « S. Lai 
University of California - Berkeley 

IBM-PC version by 
Oeotech International 
(312) 939-T16E 
May 1991 


input file name » kbnSOk . in 
Output file name * kbnSOk.ot 
Punch file name » fds 

Date and time of run •• l/Ot/1980 9:47 


MAX. NUMBER OF TERMS Ifi FOURIER TRANSFORM « 4096 
NECESSARY LENGTH OF BLANK COMMON X « 2561 9 
EARTH PRESSURE AT REST FOR SAND * .500 

OPMON 8 READ RELATION BETWEEN SOIL PROPERTIES AND STRAIN 


CURVES FOR RELATION OF STRAIN VERSUS SHEAR MODULUS AND DAMPING 


MODULUS AND DAMPir4G VALUES ARE SCALED FOR PLOTTING 


SHEAR MODULUS CLAY > HAT 1972 
DAMPING CLAY - MAY 24, I V7£ 

SHEAR MODULUS SAND SO. ROOT REL . ~ MAY 1972 
DAriPn4G SAND - FEBRUARY 1971 
ATTENUATION OF ROCK AVERAGE 


DAMPING IN ROCK AVERAGE 

9/4 




Cl ( CL A Y P I » O'- 1 0 ) M OD UL U S 
DAMPING CLAY - MAY ?.4 , 1972 

REDUCTION 
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FEB 
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MULTIPLICATION FACTOR - 
MULTIPLICATION FACTOR • 
MULTIPLICATION FACTOR » 
MULTIPLICATION FACTOR « 
MULTIPLICATION FACTOR = 
MULTIPLICATION FACTOR » 
MULTIPLICATION FACTOR « 
MULTIPLICATION FACTOR » 
MULTIPLICATION FACTOR « 
MULTIPLICATION FACTOR » 
MULTIPLICATION FACTOR « 
MULTIPLICATION FACTOR « 
MULTIPLICATION FACTOR » 
MULTIPLICATION FACTOR ■« 
MULTIPLICATION FACTOR - 
MULTIPLICATION FACTOR • 
MULTIPLICATION FACTOR « 
MULTIPLICATION FACTOR *• 
MULTIPLICATION FACTOR « 
MULTIPLICATION FACTOR « 
MULTIPLICATION FACTOR » 


. 01 OOO 
1 . 00000 
.50000 
1 . 00000 
, 01 OOO 
1 . OOOOO 
30.00000 
1 . OOOOO 
30.00000 
1 .OOOOO 
30.00000 

1 . 00000 
30.00000 

1 .00000 
30.00000 

1 . OOOOO 
30.00000 
1 . 00000 

30.00000 
1 . OOOOO 

30. 00000 
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by 

Purdue University 


•“•“Slope Stability Analysis — 
■Simplified Janbu, Simplified Bishop 
or Spencer's Method of Slices 


Run Date: 

Time of Run: 

Run By : 

Input Data Filename: kS.in 

Output Filename: k2.ot 


PROBLEM DESCRIPTION PCSTABLS Example Problem (Sarni STPP) 
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£6 

1 65 

. 00 

55. 

, 00 

£7 

1 6S 

. 00 

55 

. 00 

£S 

1 0 0 

. 00 

50 

. 00 

£9 

1 1 0 

, 00 

45 

. 00 

30 

1 1 0 

. oo 

45 

. 00 

31 

£ 0 

.00 

SO 

. 00 


X'"Right 

Y-Right 

Soil Type 

( ft ) 

(ft) 

Below Bnd 

20 ,00 

o 

o 

o 

(U 

1 

80 , 00 

50 .00 

1 

100.00 

50 .00 

1 

130,00 

60 . 00 

2 

150.00 

60 . 00 

2 

180.00 

70 .00 

2 

£O0 . 00 

70 . 00 

£ 

£30 . 00 

80 . 00 

2 

£50 . 00 

80 .00 

S 

£80 , 00 

90 .00 

2 

300 , 00 

90 . 00 

2 

330.00 

100.00 

2 

350 . 00 

100.00 

2 

"365.0 0 

■ 95.00 

2 

' 465 . 00 

95.00 

3 

395.00 

85.00 

2 

315.00 

85.00 

2 

395.00 

85.00 

3 

345.00 

*75.00 

£ 

345.00 

75 .00 

3 

£95.00 

65 . 00 

2 

£15,00 

65.00 

£ 

£95.00 

65.00 

3 • 

£45.00 

55,0 0 

2 

165.00 

55,00 

2 

£45.00 

55.00 

3 

195.00 

45,00 

2 

110.00 

45,00 

1 

195.00 

45.00 

3 

160.00 

20.00 

1 

400 , 00 

20 . 00 

1 



I SOI R OP 1C SOIL PARAMETERS 


3 Typ'H s ) of Soil 


Soil 

Total 

Sat ur at ed 

Cohesion 

Frict ion 

Pore 

Pressure 

Piez . 

T y p e 

Unit Wt 

, Unit Wt . 

Intercept 

Angle 

Pressure 

Constant 

Surface 

No , 

(pcf ) 

(pcf > 

(psf ) 

( deg ) 

Param . 

(psf) 

No . 

1 

107.0 

1 IE . 0 

. 0 

30.0 

. 00 

. 0 

1 


91.0 

96.0 

. 0 

35.5 

. 00 

, 0 

1 

3 

80 . 0 

85 . 0 

. 0 

25 . 0 

. 00 

. 0 

1 


1 PIEZOMEfRIC SURFACE(S) HAVE BEEN SPECIFIED 

Unit Weight of Water » 62.40 


Piezometric Surface No. 1 Specified by 14 Coordinate Points 


Point 

X ~W at er 

Y-Wat er 

NO . 

(ft) 

( ft ) 

1 

£ 0.0 0 

£0 .00 

£ 

3 0.00 

20 . 00 

3 

1 1 0 . 0 0 

4*5 .00 

4 

13 0 . 0 0 

45. 00 

S 

165.00 

55.00 

e 

185.00 

55 . 00 

7 

£15.00 

65. 00 

S 

2 35 . 0 0 

65 . 00 

■n 

QO 

75.00 

1 0 

2 85 . 0 0 

75.00 

1 1 

3 1 8: , 0 0 

85,00 

1 £ 

3 35 . 0 0 

85 . 00 

T 3 

305 . 00 

95. 00 

1 4 

465 , 00 

95. 00 


A Critical Failure Surface Searching Method, Using A Random 
Technique For Generating Circular Surfaces, Has Been Specified. 


too Trial Surfaces Have Been Generated. 


10 Surfaces Initiate From Each Of 10 Points 
Along The Ground Surface Between X « 10.00 

and X =. 30.00 


Equally Spaced 
ft , 
ft . 


Each Surface Tervninates Between 

and 


X » 330 .00 ft . 
X » 3S0 .00 f t . 


Unless Further !. imitations Were Imposed, The Minimum Elevation 



At Which A Surface Extends Is Y = 


.00 ft. 


8.00 ft. Lin-e S»gtT)f?nts Define Each Trial Failure Surface. 


Following Are Displayed The Ten Most Critical Of The Trial 
Failure Surfaces Examined. They Are Ordered - Most Critical 
First. 


Y’ » Safety Factors Are Calculated By The Modified Bishop Method » * 


Failure Surface Specified By 44 Coordinate Points 


Point 

X-Surf 

Y-Surf 

No . 

( ft 3 

( ft ) 

1 

23 . 33 

21 . 67 

2 

31.23 

£0 . 36 

3 

39.14 

19.20 

4 

4f . 06 

16.19 

5 

55. 03 

17.34 

0 

5 3.00 

15-64 

T 

70 , 08 

15.09 

8 

76 . 97 

15.70 

9 

66 . 97 

15.45 

1 0 

94.97 

15.35 

1 1 

102.97 

15.43 

1£ 

•M 0 . 9 r 

15.54 

13 

116.95 

16.01 

14 

126.94 

16,54 

1 5 

134.91 

17.21 

1 & 

1 4 2 . 67 

16.04 

IT 

150.61 

19,02 

IS 

1 56 . 73 

£0.15 

1 V 

150. 52 

21.44 

20 

174.49 

22.67 

21 

162 . 34 

24.45 

22 

190.15 

26.19 

£3 

1 97 . 92 

28. 08 

24 

205 . 55 

30.11 

25 

21 3 . 36 

32.29 

25 

£21 .01 

34 . 52 

2T 

228. 62 

37.09 

2S 

236. 18 

39.71 

20 

243 . 68 

42. 47 

30 

c55 1 , 1 4 

45,38 

31 

256.53 

48.43 

32 

255.67 

51 .61 

33 

273. 15 

54. 94 

34 

280.35 

58.41 

35 

267.50 

52. 02 

35 

£94.57 

55.75 

37 

301.57 

59. 53 

38 

306 . 49 

73. 54 

3? 

3.15. 33 

TT. 78 

40 

322. 1 0 

82. 05 

41 

326 . 78 

85 . 45 

42 

335 . 37 

90. 98 

4 3 

341. 68 

95. 54 

41 

34 7 75 

100.00 



C 1 r- c I e C e n 1 e r At X 


9S.6 ; Y - 433.2 and Radius 


41T.8 


♦ ^ t . S62 ^ 


Failure Surface Specified By 44 'Coordinate Points 


Point X-Surf Y-Surf 

No . ( ft ) Cft ) 


! 

23. 

33 

21 . 

67 

2 

31 . 

20 

20. 

22 

3 

39 . 

1 0 

15. 

92 

4 

4t . 

02 

1 7 . 

79 

S 

54 . 

96 

1 6 . 

81 

0 

62. 

91 

15. 

99 

T 

TO . 

89 

15. 

34 

8 

78. 

87 

1 4 . 

84 

9 

2 6 . 

3 6 

1 4. 

50 

1 0 

■ 9 4 . 

86 

1 4. 

32 

] 1 

102. 

86 

1 4 . 

30 

T ,2 

no. 

8 6 

1 4 . 

44 

1 3 

118. 

85 

14. 

75 

1 4 

1 26 . 

84 

1 5 . 

21 

1 S 

134. 

82 

1 5 . 

83 

1 0 

MS . 

75 

16. 

61 

1 T 

15 0. 

72 

17. 

55 

18 

1 5 8 . 

65 

15, 

65 

1 9 

166 

55 

t 9 , 

. 90 

2 0 

1 T4 . 

<12 

21 . 

. 32 

21 

1 82 . 

, 27 

22 : 

. 89 

22 

190. 

, 08 

24, 

. 62 

2 3 

1 97 , 

. 85 

26 

. 50 

24 

205 

, 59 

28 

. 54 

2S 

2 13 , 

. 28 

30 

. 73 

26 

220 

. 93 

33 

. 08 

2T 

VfWp 
W V 

. 53 

35 

. 58 

t? o. 

236 

. 08 

38 

. £3 

29 

24 3 

. 57 

41 

. 03 

30 

25 1 

.01 

43 

. 98 

31 

258 

. 38 

47 

. 08 

32 

265 

. 69 

50 

. 33 

3 3 

272 

. 94 

53 

. 72 

3 4 

28 0 

. 12 

57 

. 25 

35 

287 

. 22 

60 

. 93 

36 

29 4 

. 25 

64 

. 75 

37 

301 

. 20 

68 

.71 

38 

308 

. 07 

72 

. 81 

39 

31 4 

. 36 

77 

. 04 

40 

321 

. 56 

SI 

.41 

41 

328 

. 17 

S5 

. 92 

4 2 

334 

. 69 

90 

. 55 

43 

341 

. 12 

95 

, 31 

44 

34 7 

19 

1 00 

. 00 


Circle Center At X = 9^,8 j: T « 414.2 and Radius^ 


399,9 


,T 


1 .S(:,S 
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